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Currents of the single-electron transistors driven by time-dependent fields via external dissipative 
circuits are investigated theoretically. By expressing the external circuit in terms of driven har- 
monic oscillators and using the reduced-density operator method, we derive time- and environment- 
dependent tunneling rates in the regime of sequential tunneling and present expressions for both 
displacement and tunneling currents with these tunneling rates. It is found that the dissipative 
environments affect tunneling currents in two ways; the determination of driving voltages at tun- 
neling junctions and the depletion of particle-hole distribution functions. Considering a simple 
dissipative circuit, we discuss the effects of the environment on tunneling currents in both static 
and time-dependent cases. 
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I. INTRODUCTION 

There have been considerable interests in a single- 
electron transistor because of its potential applications. 
Since Coulomb blockade peaks in its conductance oscil- 
lation are very sensitive to a fraction of charges, it is ex- 
pected to be one of promising candidates for a detector 
in measuring quantum states of quantum computation 
and information processing |lj. Recently, Schoelkopf et. 
al.^ developed a high sensitive single-electron transis- 
tor. By introducing a LC-resonant circuit connected to 
a single-electron transistor, they drove the system in the 
radio-frequency regime to overcome 1// noise and obtain 
a high sensitity of detecting charges. 

From the theoretical point of view, such a single- 
electron transistor is also very interesting because 
one should consider the influence of dissipative 
environments (LC-resonant circuit with a cable resis- 
tance) on current-voltage characteristics as well as effects 
of time-dependent external perturbations. In the pres- 
ence of the dissipative environments, tunneling rates of 
quasiparticles are strongly affected because an additional 
energy is needed to excite the environments as shown in 
the case of single- and multi-junctionsQ. On the other 
hand, time-dependent perturbations affect phase coher- 
ence of quasiparticles in time and its effects manifest itself 
in photon-assisted tunneling; quasiparticles are also able 
to tunnel by absorbing or emitting photons0,|H. Thus, 
quasiparticles in the single-electron transistor have two 
different energy-exchange mechanisms when they tunnel 
through junctions and then, one should take into account 
these two mechanisms in calculating the tunneling rates 
simultaneously. 

In this work we investigate effects of the dissipative 
environments and the time-dependent perturbations on 
current-voltage characteristics of a single-electron tran- 
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sistor by calculating the tunneling rate and present ex- 
pressions for currents at each electrode. In calculating 
the tunneling rate we make two assumptions for our 
model of the single-electron transistor. Firstly, tunnel- 
ing barriers are considered so opaque that quasiparticles 
in each electrode are well localized there and their mo- 
tions can be described by separated Hamiltonians. In 
other word, the tunneling barriers have a resistance much 
larger than the resistance quantum Rk — h/e^. Sec- 
ondly, it is assumed that time between successive tunnel- 
ing events is much larger than the charge relaxation time 
of the dissipative environments. This assumption makes 
the problem easy to treat the dissipative environments as 
heat reservoirs being in thermal equilibrium. However, 
we still treat quasiparticles confined in the region of a 
quantum dot as in non-equilibrium. 

The paper is organized as follows. We first describe the 
Hamiltonian of the single-electron transistor in Section 
II. By expressing the dissipative environments in terms 
of driven harmonic oscillators, we give the Hamiltonian 
separating it into an interested system part and its en- 
vironment, and derive time- and environment-dependent 
tunneling rates using a reduced-density operator method. 
In Section HI, we present expressions for currents flowing 
in each electrode in terms of displacement and tunneling 
components. As an applilication of our expressions, tun- 
neling currents driven by a simple dissipative circuit are 
examined in Section IV, and then, a brief summary is 
given in Section V. 



II. HAMILTONIAN 

In Fig. 1, we show a typical drawing of the single 
electron transistor driven by time-dependent voltages via 
possible dissipative elements. The Hamiltonian of the en- 
tire system is modeled by H = 'Hqp+'HRLc{t)+'HT where 
the first two terms describe the motion of the system in 
the absence of tunnelings and the last one is the tunnel- 
ing Hamiltonian. In the absence of tunneling, the system 
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FIG. 1: A typical drawing of the single-electron transistor is 
shown where time-dependent voltages are applied to a quan- 
tum dot via possible dissipative elements connected to drain 
and gate electrodes, respectively. 



is viewed as just a simple electronic circuit(a lumped- 
circuit) because all tunneling junctions are considered as 
capacitors. Then, its motion is described by two inde- 
pendent degrees of freedom; microscopic and macroscopic 
variables 0. 

By the microscopic variables, we mean those of quasi- 
particles in conductors of Fig. 1. In our system, since 
quasiparticles relevant to tunneling are distributed in 
a island called a quantum dot and adjacent electrodes 
(source and drain), we denote their Hamiltonian as 



qp 



'Hdo 



Hiead where Hdot = J2i ^idj^i + Vee{dl,di) 
and Hiead = J2a=s..DJ2k^kaal^aka- Here, di and aka 
{dj and a^^) are the annihilation(creation) operators for 
quasiparticles in the quantum dot with an energy e; and 
in the electrodes a with an energy eka, respectively. In 
the quantum dot we denote the electron-electron interac- 
tion by Vee{dl , di) which in the most case is approximated 
by a simple Coulomb-blockade model Q^/2C^ with ex- 
cess charges Q and — Cd + Cs + Cg- In the metaUic 
electrodes, energies of quasiparticles are assumed to be 
measured relative to their chemical potentials and in- 
dependent of time-dependent external perturbations by 
considering small frequencies not to excite plasmon. 

On the other hand, macroscopic variables represent 
charges on capacitors and flux through inductors in the 
lumped-circuit. To describe their motion, we model the 
dissipative elements of Zo{uj) and Zcii^) in Fig. 1 with 
different sets of simple harmonic oscillators {i.e. LC- 
circuits) following Caldeira and Leggettj]^. Then, start- 
ing from the Lagrangian formulation where generalized 
coordinates are chosen to be 0(i) = J dtev{t) /Hior a, volt- 
age v{t) across each capacitor, the Hamiltonian TlRLc{t) 
of the macroscopic part is expressed as. 
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where charges Qj and (?„ are conjugate to phases h4>j/e 
and hipj/e, respectively, and satisfy the commutation re- 
lations of [4>j,Qj] = ie and [ipn^qn] = ie. The external 
perturbation Va (t) are treated as a classical field and in- 
volved in externally driving phases ipoit) and ipcit) as. 
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h Cd 
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(2) 



The first two terms in Eq. are additional charging 
energies on the capacitors connected to each electrode 
with effective capacitances Ci = CdCs/{Cd + Cs) and 
C2 = C1,Cg/C-e{Cd + Cs), whereas the remaining terms 
describe charging energies on capacitors with a capaci- 
tance Cm and magnetic energies of inductors with a induc- 
tance Lm in the dissipative elements, Z-oiuS) and Zq(ijj). 
In fact, the charges Qi, Q2, and Q in Hrlc describe 
charges Qa (a = S, D, G) accumulated on a plate of the 
capacitor connected to the electrode a in Fig. 1 and they 
are related to each other by. 
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Then, classical relaxation equations for the charges Qa {t) 
are exactly recovered from the Heisenberg equation of 
motions. To do this, the impedances Zq(uj) and Zpiuj) 
at a angular frequency lo are anticipated to have sets of 
parameters {cn,Ln} and {cm,Lm} fikeQ, 
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The tunneling part of the Hamiltonian Tif may be 
given by. 



akl 



(5) 



where denotes the element of tunneling matrix be- 
tween a state | k) in the electrode a and a single particle 
state I I) in the quantum dot. Here, the phase operator 
4>a (ck — S, D) is defined in terms of <f>i and (/)2 through 
the relation. 



= A 



4>2 




(6) 



Thereby, the operator e"''° increases excess charges by 
the elementary charge e in the electrode a for every tun- 
neling event because they satisfy the commutation rela- 
tion [0Q, Qq] = ie. 
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In reality, the form of the Hamiltonian 7i = Ti-qp + 
T^RLc{i)-^'^T is not adequate in obtaining expressions of 
currents because the time-dependence of Ti-RLcit) makes 
the problem complicated. To circumvent this, we per- 
form a standard time-dependent unitary transformation 
with, 

2 N+M 

n n exp{ixj{t)(j).j +ix'^{t)ip,n} 

j — l m— 1 

exp{-mj(t)(3j/e - iu'^{t)qm/e}. (7) 

Here, the sets of time-dependent functions {xj{t), x'j{t)} 
and {umit) , u'^{t)} are chosen to rotate the Hamiltonian 
T^RLcit) into the form irrelevant of time, TLRLci^)- By a 
straightforward calculation for the rotated Hamiltonian 



■HR{t) = U\t)HU{t) - im^'it) 



dUjt) 



(8) 



one can show that a time- independent form of Ti.RLc{0) 
is obtained by choosing Uj{t) just as the phase difference 
across capacitors forced by ipoit) and ipcit) in the ab- 
sence of tunneling and Xj{t) = Cjhuj/e^. That is, Uj{t) 
is chosen by, 



ui{t) 

U2{t) 



= A 



VD{t) 
Vs{t) 



I dr[{5^^n -^)V'r>-^VB + ff.(r)} (9) 



where we separate the external perturbation into DC and 
AC parts, i.e., Va[t) — V^ + Va{t). Here, ga[T) describes 
voltage differences at tunneling junctions exerted by the 
AC part of the voltage sources and is given in its Fourier 
component as. 




vd/Zd 
vg/Zg 



(10) 



where Va{(jj) is a Fourier component of the AC part of 
the external perturbation and the matrix Z is defined as, 



Z = A 



iuiCi + Zj^ 



ZLjCz + Z'r 



A^(ll) 



Finally, under the unitary transformation U{i), the total 
Hamiltonian becomes, 

Wfl(t) = 7^rfot+H/ead+H7?Lc(0)+H^(t). (12) 

Here, H'rp{t) is the tunneling Hamiltonian rotated by the 
unitary transformation U and now involves effects of the 
time-dependence perturbations as. 
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A. Ensemble average 

In order to evaluate an ensemble average of phys- 
ical quantities, we use the reduced-density operator 
method^. In this method, the entire system is divided 
into an interested system part being in non-equilibrium, a 
heat bath, and their interaction; Hr = Tis +'Hb +'Hsb- 
Then, the effective density operator describing the in- 
terested system is obtained by averaging the Liouville 
equation over the heat bath. Since the heat bath is con- 
sidered as in thermal equilibrium, the ensemble average 
of physical quantities is expressed as the sum of expec- 
tation values between quantum states weighted with the 
reduced-density matrix. 

In our case, we consider the system part as quasiparti- 
cles in the quantum dot, Tis = Ti-dot, and the heat bath 
as the combination of quasiparticles in the electrodes and 
the lumped-circuit, i.e. Hb = 'Hiead + 'Hrlc{^)- Then, 
the system and the heat bath are coupled by the time- 
dependent tunneling Hamiltonian HsB{t) — 'H'rp{i). By 
adopting very opaque tunneling junctions (weak coupling 
of Ti-SB) such that a typical tunneling period is much 
larger than the charge relaxation time of the heat bath, 
one can treat the heat bath as being in thermal equilib- 
rium. In this case, its density operator is proportional 
to e~^'^^ at an inverse temperature (3. 

The density operator describing the interested system 
p{t) is obtained by tracing the density operator ptot{t) for 
the entire system over the heat bath, p[t) = tT:B{ptot{t)}i 
and its equation of motions is derived from the Liouville 
equation. In a non-Markovian form, p{t) is given by0, 



(14) 



Here, the generalized scattering operator C{t) stands for, 
= h j dTtrB[[nsB{t),W(T,t)pB,H'sB{r,t)'^] 



-o{n 



SB) 



(15) 



where p^ (t) and T-Cgg{t) are their interaction pictures of 
p{t) and TCsB{t), respectively; 



p\r,t) = e 
n'sBir,t) = e 



p{T)e 

-i{ns+HB){t-T)/h 

J{Hs+HB){t-T)/h 



ms(t-T)/h 



nsBir) 



(16) 



(17) 



Here, the forth order contributions of the interaction 
Hamiltonian TLsB{t) are ignored by taking into account 
opaque tunneling junctions. Then, with this density op- 
erator, one can express the ensemble average of an arbi- 
trary operator O in terms of relevant system operators 
by replacing the part of the heat bath with their equi- 
librium values. That is, (O) — trs{trs{Optot(i)}} = 
trs{Osp{t)} where trg means the average of the opera- 
tors over the system and Os is a relevant system operator 
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B. Expectation values for heat- bath operators 

For the ensemble averages of physical quantities, 
one need to evaluate several time-correlations between 
heat-bath operators. First, for those of the RLC- 
circuit, the correlation between phase operators of 
trB{e~"^°e*'^°(*)ps} is necessary to calculate its effect 
on tunneling. Here, (t>a{t) is the Heisenberg operator of 
(j)a with respect to TLrlc- Since the iiLC-circuit of Eq. 

is considered as the sum of independent harmonic os- 
cillators in equilibrium, this correlation function can be 
rewritten as[ 



trB{e 



Pb} 



,trB{(0„0„(t)-0„0„)pB} 



(18) 



Then, based on the linear response theory, one can 
show that the fluctuation of Ja{t) = tTB{{4'a4'a:{t) — 
4>a4'a)pB} is directly related to the dissipation of the 
-RLC-circuit (Fluctuation-dissipation theorem). From re- 
sponse functions of the i?LC-circuit of Eq. iQJ together 
with Eqs. © and lO, Ja{t) is given by. 



Lo Rk 1 - e-^!^ 



(19) 



where Z^{u;) is an effective impedance of the i?LC-circuit 
seen from the tunnel junction and is equal to diag- 
onal elements of the impedance matrix Z of Eq. (|ll|l : 
Z\~i{uj) — Zii{uj) and Zg(uj) = Z22(w), respectively. 
In reality, Zj^(uj) has a slightly different form from the 
impedance seen from a tunneling junction a in Fig. 1. 
This is because the region of the quantum dot in Fig. 1 is 
independent of the i?LC-circuit and does not contribute 
to the fluctuation in the absence of tunneling. 

As for the electrodes, the following particle and hole 
evolutions are necesarry to evaluate tunneling currents, 

ti-B{ai^akc{t)pB} = /FD(efea)e^^"l*l-^^'="* 
trsiflfcaaLW-os} = {l-/FD(efca)}e-^°l*l + ^'''°(^0) 

where a^ait) is the Heisenberg operator of a^a with re- 
spect to Hiead and fpoif) = 1/(1 + e''') is the Fermi- 
Dirac distribution function. Here, the exponential decay 
of 6"'''° 1*1 represents effects of tunneling on states in the 
electrodes. In reality, the heat bath as well as the quan- 
tum dot in our system are affected by tunneling. As a 
result, the evolutions of quasiparticles in the heat-bath 
are different from those in the isolated one and are usu- 
ally represented by a finite life time. The exponential 
decay is inserted by hand to represent such a effect on 
states 0. The parameter "fa is assumed to be the bare 
tunneling rate at the electrode a, 



I Tt, p 5{e - eua), 



(21) 



which is usually considered as a constant independently 
of an energy within a so-called wide-band limit; ia{e) = 

la- 



C. Time- and environment-dependent master 
equation 

Now, we evaluate the reduced-density operator of Eq. 
(|14|l in the basis representation. A simplified form is 
obtained when one expands the reduced-density operator 
in terms of many-body eigenstates(| r) and \ s)) of the 
quantum dot, i.e.. 



p(i)=^P..(Ok}(^ 



(22) 



Substituting this into Eq. (|14|l and then projecting it on 
one of diagonal components, the occupation probability 
Pss{t) at a state | s) is given as a balanced form. 



dt 



ra^— +, 


dTP,,(t + T)r^«(t,T) 



(23) 



where the first term describes the increasing rate of the 
probability density by transitions from other states while 
the second term is a decay rate due to transitions to oth- 
ers. In deriving this result we disregard the contribution 
of off-diagonal components because their effects are the 
forth order of the interaction Hamiltonian, 0{T-Cg^) and 
thus, the result does not contain coherent evolutions be- 
tween many-body states which may be caused by external 
perturbations. 

In Eq. H23() the memory kernels of r"^^(t, r) describe 
quasiparticle tunneling into (-I-) or from (— ) the quantum 
dot through the barrier a to result in the transition from 
a state | r) to another state | s). The detailed forms of 
r'^^{t,T) are given by. 

Re r°° 

KHt,r) = — deA'^^{e + E^,^) 

{%€ ZG "1 

laT -j-'r + J„(±t) - - dt'ga{t')j (24) 

Here, A'^^{E) represent tunneling rates of quasiparticles 
into(+) or from(— ) the quantum dot with an energy gain 
of E when there is neither dissipative elements nor alter- 
nating perturbations in addition to negligible collision- 
broadening. In this case, the tunneling rates A"^{E) are 
reduce to the widely used formulajl^ ; 

A'^HKs^) - 7o5± {i T ^ ± fFDiE^s^)} (25) 

where coefficients 5^ represent the selection rules of tun- 
neling, 

S+ = 5] |(r I I \\ = ^ |(. I di I r) 1^(26) 

and = ±{E, - Er) + - %^)eV^ - %ieVS, 

are energy gains at a tunneling event with Er and _Es, 



5 



eigenenergies of states | r) and | s), respectively. On the 
other hand, the exponential part of the integrand in Eq. 
H24|l includes the effects of the dissipative environment 
and the alternating perturbations. As indicated by the 
Ticn-Gordcn theory 11], the external source contributes 
the imaginary part in an argument of the exponential 
function via ga{t), and thus affects phases of electronic 
states. Whereas, the term Ja{t) in general have both real 
and imaginary parts, and thus gives rise to the damping 
of states as well as the change of phases. Roles of Ja{t) 
in r"^^ {t, t) are more easily understood by transforming 
into its Fourier components as, 



(27) 



Substituting this relation into Eq. H24|l . one can see that 
tunneling occurs with a weight of Pa{i^) in the range 
between uj and uj + dto and corresponding energy gains 
become e ± hoj + E"^. Since the sum of Pa{uj) over all 
frequency range is equal to one, Pq(cj) can be interpreted 
as the probability density to exchange the energy huj be- 
tween the system and its environment 6j|. By considering 
energy differences between tunneling events, one can see 
that Paii^) in the positive (negative) frequency represents 
the probability to emit(asorbe) photons to(from) the en- 
vironment. 

From Eq. (|24|l . it is noted that effects of the dissipa- 
tive environments on tunneling are two folds. The first 
is the determination of the voltage difference gait) at 
each tunneling junction through Eq. (|10|l . The other is 
the probability density Pa{ui) which is determined by the 
characteristic impedance (w) of the environments via 
Eq. (O. 



III. EXPRESSIONS FOR CURRENTS 

Now, we calculate currents in the electrodes in Fig. 1 
where a positive current at each electrode is defined to 
flow into the quantum dot. The current flowing in each 
electrode a consists of two different contributions; tun- 
neling currents of quasiparticles (t) and time- variation 
of charges on capacitors in the lumped-circuit /^(t) called 
displacement currents, i.e., 



ut) = ii{t)+iiit). 



(28) 



Here, each component is calculated by time-derivatives 
for the ensemble average of particle numbers and charges; 



dt 



dt 



{Qa)o 



(29) 



where (. . .)o means the average over the Hamiltonian of 
TC. With the above current expressions, the currents flow- 
ing into the quantum dot are conserved even for the sys- 
tem subject to the time-dependent perturbations, as em- 
phasized by Biittiker in his recent work|3|. This can be 



FIG. 2: For calculating both displacement and tunneling cur- 
rents in the single-electron transistor of Fig. 1, an equivalent 
circuit is shown, where components of tunneling currents are 
considered as current sources. 



shown by calculating the time-derivative of (Q), which 
gives d{Q)/dt = {[Q,'H\)/ih = 4(t) + Ijjit) reflecting 
the fact that the increase of charges in the quantum dot 
is enabled by tunneling processes. Alternatively, since 
d{Q) /dt is the sum of all displacement currents out of 
the quantum dot from Eq. (|3Jl, we obtain the conserva- 
tion of the currents, /^(i) = EMit) + lUt)} ^ 
with /^(i) = 0. 

By solving the Heisenberg equations of motion, 
d{Qa)/dt = {[QcHD/ih, we can express the displace- 
ment currents in terms of the contributions of tunneling 
currents and external perturbations. The results are, in 
its Fourier components of /^(w), 



1 



Zg 



iuCg 

vd (w) 
-vg{uj) 



iujCs 



Z 



D 







iluCg 



iuCg 



where lai^) is a Fourier component of a tunneling cur- 
rent Ia{t). Here, the first term on the right-hand side 
is the contribution from alternating perturbations while 
the second terms are resulted from tunneling. Alterna- 
tively, the above results can be expressed in the equiva- 
lent circuit as shown in Fig. 2 by modeling the tunneling 
contributions as current sources. Then, once tunneling 
currents /^(i) are known, the total currents at each elec- 
trode are determined by applying basic circuit rules to 
Fig. 2. 

Now we evaluate the tunneling component of the to- 
tal current in Eq. I|28(l . From the Heisenberg equation 
of motion, d{eJ\fa)o/dt = ([eA/'a, W]}o/i^, the timneling 
current /^(t) is calculated as. 



H.c. 



akl 

akl 
(Ia(i)) 



H.C. 



(31) 



where in the second line (...) denotes an ensemble av- 
erage over the rotated system Tin. By considering the 
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second order contribution of 7is b or only sequential tun- 
neling processes, the expectation value of Eq. can 
be rewritten in the reduced-density operator formalism 



as, 



Here, comparing this with Eq. one can see that 

the part of the commutator is resulted from the evolu- 
tion of the total density ptot{t)- By replacing heat-bath 
operators with their equilibrium values and substituting 
the reduced-density operator in the basis representation, 
finally we obtain the tunneling current at the electrode 
a as, 

As in the case of the occupation probabilities, off- 
diagonal contributions are also ignored because their ef- 
fects are the forth order of TLsb- 

Equations ^'2.'A\ . 1)24(1 . and H33I) are main results of our 
work. Based on the master equations of Eq. (|23|l and 
the memory kernel of Eq. H24|) . the displacement and 
tunneling currents flowing in each electrode of Fig. 1 can 
be calculated using Eqs. H3U|I and (|33|l . respectively. 



A. Expressions in the Fourier space 

For numerical calculations, it is useful to introduce the 
Fourier transform of the memory kernel, Eq. (|24|l . When 
external perturbations are periodic functions with an an- 
gular frequency uja, we define Fourier components as. 



It: I UJA 



dt 



dr 



exp { — ilujAt + imujAT}T^s{t,T). (34) 

Then, by writing the external perturbations in their 
Fourier components, 



exp{ 



dt'g^it')}^ J2 e""^*77aW, (35) 



n— — oo 



we obtain ^^^{l.m) as, 



(36) 



Here, effective quasiparticle(+) and hole(— ) distributions 

hr. 



f^^{E) are defined by. 



dljjPa{uj) 



27ri 



-I, —my 



(37) 



where 



1-K 



2ni 



{E — (m + n)hujA i hio} 



and i}^{z) is a digamma function. Without the dissipa- 
tive elements and the alternating perturbations in ad- 
dition to negligible collision-broadening 7^, f^^{E) are 
reduced to ^{1 T tanh(^)}, the Fermi-Dirac distribu- 
tions for particles and holes(— ). However, in general, 
fim deviated from the Fermi-Dirac form due to the 
dissipative environments(described by Pa) as well as the 
finite life time of quasiparticles (described by 7^). 

In terms of r^J^(/, m), the occupation probabilities and 
tunneling currents are calculated as, 



iIujaPss{1)^Y. J2 {Prrim)t^l{l~m,m) 

5(m)f^«(;-m,m)} (38) 



and 



(0 = e E Prr (m) {f^?^ (/ - m, m) - C+ {l-m,m) )(39) 

rsm 

where we also expand Pgg (t) and (t) in their Fourier se- 
ries; P,s{t) = E; e*'"-"*Ps.(0 and /* (i) = e*'"^*J* (0- 
In special case of Pai^^) = S{uj), i.e. without dissipative 
elements in the circuit, we find that the above results to- 
gether with the tunneling rates of Eq. (|35|l are similar to 
the formalism developed by Bruder and Schoeller,^. 



B. Time-convolutionless form 

Evaluating the integral of Eq. 124(1 , we derive another 
form of the memory kernel, 



csch(f5) 



r?.^(i,r)=7„5± 5{t)t 



hp 



lE^^T/h+ Ja{±T) 



dt'gait')} (40) 



t + T 



with r < 0. Since this function decays exponentially from 
T = 0, we now expand Prr(t + T) in Taylor series at t = 
to calculate Eqs. 1(23(1 and 1(331) . By collecting the leading 
contributions, the occupation probabilities and tunneling 
currents are shown to be. 



dPssjt) 
dt 



J2 [Prr{t)KHt) ^ P,,{t)T^Ht) 



+o(h|b) 



(41) 



and 



m = eY,P,rit){T^rit)~K+it)} (42) 
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where T^l{t) = drF^Ki, r). The next contribu- 
tions are the fourth order of the interaction Hamilto- 
nian Ti.sB{t) (These results can be also derived start- 
ing from the time-convolutionless solution of the density 
operator 9J) and, neglecting them in the spirit of sequen- 
tial tunneling, the results are now time-convolutionless. 
In the Fourier space, the time-convolutionless results read 
as, 

ram ^— -|-.— 

-A,(m)f^«(Z-m,0)} (43) 

and 

rsm 




FIG. 3: We show a simple example of the single-electron 
transistor driven by a time-dependent field via a resistor Rq 
in the drain. Here, symmetric tunneling barriers are assumed 
with their bare tunneling rates of 70 = 7s = 7o/2 and capac- 
itances of Cs = Cd = 10~*e^/fi7Q. The charging energy is 
chosen to be Ec = 100^70 and all numerical calculations are 
done at a temperature of fcsT = ^70. 



C. Adiabatic limit 

As shown in Eq. H40(l . r"^{t,T) is dominant around 
T = and thus, for a slowly varying external field {uja ^ 
7a + TT /hp) it can be further approximated as. 



r.T(i,T) 



<5(r)T 



CSCh{7TT / hp) 

h(3 



f *e 

exp <^ 7ar - — r + J, 



a(i)- 

a(±r)} 



(45) 



Namely, in this limit, energy states in the quantum dot 
are merely modulated by the external perturbations. 
Furthermore, since these modulations are much slower 
than the equilibrated rate of 7q + n/hp due to tunneling 
and temperature, one can treat the problem as a static 
one with an additional DC-bias voltage of ga{t) at each 
instant. Then, the occupation probabilities may be de- 
termined from a static balance relation, 



0= Yl [Prr{t)KHt) - Pss{t)T'^Ht) 



(46) 



IV. APPLICATIONS OF FORMALISM 

As applications of our results, we now examine time- 
dependent currents in a single-electron transistor based 
on Eqs. and To understand I - V 

characteristics easily, we consider a simple circuit of Fig. 
3 where a sinusoidal voltage of voit) — v^cos{uj At) is 
driven at a drain through a resistor Rq and a static volt- 
age is applied on a gate without any impedance. The 
electron-electron interaction is assumed by a Coulomb 
blockade model; Vee = Q^/2Cs. For well-separated 



peaks of the Coulomb blockade oscillation, a charging 
energy {Ec = e^/2Cs = IOO/170) is chosen to be much 
larger than a broadening due to tunneling hjo and a ther- 
mal energy ksT = h-jg considered here (70 is a unit for 
a tunneling rate). We also assume symmetric barriers in 
the drain and sources, and denote their tunneling rates 
with = 7s = 7o/2. Although we assume the sym- 
metric tunneling barriers, the voltage differences ga(t) 
across each tunneling junction as well as the character- 
istic impedances Z\JyuS) are different from each other be- 
cause of a special geometry of our circuit, so called, a 
common-source geometry. In this geometry, the voltage 
difference of the tunnel junction to the source is mainly 
determined by the gate voltage Vq while the junction to 
the drain depends largely on the drain voltage Vo- 



A. Effects of the dissipative element Jio 

As mentioned in the previous section, the first role of 
the dissipative elements is the determination of the AC 
voltages (?a(i) across each tunnel junction. According to 
Eq. H1U|) , (?Q it) for the circuit of Fig 3 are given by, in 
their Fourier components. 



Cs + Cg 
-Cd 

1 

Cp 
Cg 



Vd (i^) 



Ce -I- iwiCc 



1 + iwCdRq 



Cs)CdRo 
(47) 



where a relatively large capacitance of the gate capaci- 
tor compared with those of the source and drain is as- 
sumed in the second relation. It is noted that for a small 
frequency compared to a cut-off frequency 1/CdRo, the 
most part oi vd is applied on the drain junction with a 
small fraction of Cd/Cq on the source junction. How- 
ever, above the cut-off frequency, (fa fall like 1/uj, there- 
fore, a smaller fraction of vd is applied on the tunneling 
junctions. 
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The second role of the resistance Rq appears in the 
broadening of quasiparticle distributions as in Eq. 
through the probabihty density Pa- From Eq. 
Zl^{ijj) becomes for the circuit of Fig 3, 

The impedance Zg(u}) is smaUer than Z*j^(uj) by a fac- 
tor of C|)/Cq and, for Cq ^ Cd, we can readily set 
Ps(a;) ~ 5[lo) which implies that tunneling through the 
barrier to the source is irrelevant to the environment. 
In Fig. 4, we show the probability density Pd(w) at 
the drain and corresponding particle-distribution func- 
tion for various Rq when there are no alternating per- 
turbations. Starting from a (5— function for Rq = 0, the 
probability density Pd{^) shows a Lorentzian shape in 
the region of positive frequencies and exhibit exponen- 
tially decaying behavior of Pd(— w) = exp{— fiw/3}PD(w) 
in the negative region[^. As the values of Rq increases 
the shapes of P£i{uj) are found to become more broad 
together with shifted peak positions to a positive fre- 
quency. This means that when particles tunnel a barrier 
more energies should be transferred to the environment 
as Rq increases. These results are also reflected in the 
particle(hole)-distribution hmctions as shown in Fig. 4- 
(b). As i?o increases, the particle-distribution functions 
are largely depleted in the region of negative energies. 
These depletions are similar to the case as if it is a high 
temperature and a chemical potential is shifted to a neg- 
ative energy. Thus, one can expect that, compared to 
results of i?o = 0; tunneling currents are smeared out 
and start to flow at a higher drain voltage as the resis- 
tance i?o increases. 

B. Tunneling currents under time-dependent fields 

As a simple example, we consider two degenerate 
states(ei =£2 = 0) in the quantum dot and, there- 
fore, two possible many-body states with energies of 
Eg = s^Ecis = 1,2) are available for tunneling. In the 
static case, we obtain typical I — V curves; step-like be- 
havior and Coulomb blockade oscillations as a function 
of drain and gate voltages, respectively, for _Ro = as 
shown in Fig. 5 (dashed lines). Under the dissipative 
environments (i?o 7^ 0), the shape of the I — V curves are 
found to be smeared out due to the depletion of the par- 
ticle(hole) distribution functions(not shown in Fig. 5). 

The smearing of the tunneling currents is also found 
for the case of time-dependent fields. Applying an alter- 
nating perturbation, fine structures in tunneling currents 
are developed because energy levels in the quantum dot 
are split into photon-side bands, Eg + khujA (k =integer). 
These split energy levels are well identified in a DC part 
of the tunneling current as a function of additional static 
drain voltages, especially in the case of Rq = O(dotted 
line in Fig. 5- (a)). Namely, each step appears at the 
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A (a) - 
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FIG. 4: We plot the probability density of Pd(lj) at the 
drain in (a) and corresponding particle-distribution function 
foo'^iE) in (b) for Ro = O(dotted lines), 1, 2, 3, 4, and 
5kf2(solid lines) when there are no alternating perturbations. 



voltage which gives the chemical potential of the drain 
corresponding to one of energy levels, Eg + khujA- Un- 
der the dissipative environments, this step-like behavior 
is washed out as well as the height of the steps is re- 
duced as Rq increases as shown in Fig. 5-(a)(solid fines). 
Different from the static case, the reduction of the tun- 
neling currents can be caused by the decrease of voltage 
differences across tunneling barriers as well as the deple- 
tion of the particle distributions described by Eqs. H47|) 
and H48|l . respectively. However, the results of Fig. 5 
are mostly responsible for the latter effect of the envi- 
ronment in this numerical simulation because the cut-off 
frequency of I/RqCd is still much larger than the applied 
frequency to a even for Rq = 5kf2. Thus, we find that the 
tunneling currents are smeared out nearly by the same 
amount in both static and time-dependent cases as Rq 
increases. 

For various values of Rq, we also examine a DC part 
of the tunneling current as a function of a static gate 
voltage Vq with a static drain voltage being zero, and 
plot the results in Fig. 5-(b). Instead of a Coulomb 
blockade peak in the static case(dashed line), it is found 
that the tunneling currents have negative or positive val- 
ues depending on the gate voltages, and the direction is 
abruptly altered around the peak position. This behavior 
is retained for lower frequencies of the alternating pertur- 
bation, even in adiabatic limit. In reality, the direction 
of the tunneling currents is easily inferred because the 
voltage difference across the tunneling junction to the 
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FIG. 5: When an alternating field is applied with v% — 
50/170 /e and lua = 3O70, the DC part of the tunneling cur- 
rents(in units of 670) at the drain are plotted as a function of 
a static drain voltage in (a) with Vq = 150ft7o/e and as a 
function of a static gate voltage Vq in (b) with = 0. Here, 
the thicker solid line describes the result for a larger resistance 
of Ro = 1, 2, 3, 4, and 5kf2, respectively, with the dotted line 
for i?o = 0. The step designated by kuJA is resulted from a 
photon-side band Ei + khujA- For a comparison, the results 
without the alternating field are shown in both graphs with 
dashed lines for Rq = Q where = 2/i7o/e is used in (b). 



to the source, and vice versa. Similar to the results as 
a function of a static drain voltage in Fig 5- (a), photon- 
side bands are manifested itself in steps apart from each 
other by hwA- We find that these steps are well resolved 
in the case of i?o = while they are relatively washed 
out as i?o increases. 

For valid applications of the time-convolutionless for- 
malism we also calculate tunneling currents for the circuit 
of Fig. 3 based on Eqs. (|41|l and H42|l . By varying pa- 
rameters within small tunneling rates, it is found that the 
time-convolutionless formalism give negligible differences 
from results obtained by the time-convolution forms. 
V. SUMMARY 

In summary, we have studied time-dependent currents 
of the single-electron transistor embedded in possible dis- 
sipative circuits and driven by time-dependent perturba- 
tions. In the regime of sequential tunneling, we present 
numerically tractable forms for both displacement and 
tunneling currents where the tunneling rates of Eq. I|24|) 
contain explicitly the influence of the dissipative environ- 
ments and time-dependent perturbations. We find that 
the dissipative environments affect tunneling currents of 
the single-electron transistor in two ways; the determi- 
nation of driving voltages at tunneling junctions and the 
depletion of particle-hole distribution functions at each 
electrode. Applying our formalism to a simple dissipa- 
tive system and solving the problem numerically, we show 
how steps in tunneling currents developed by photon-side 
bands are smeared out as the system becomes more dis- 
sipative. 



source is mainly determined by the gate voltage as men- 
tioned in the previous section. In the region of the nega- 
tive(positive) tunneling currents, the chemical potential 
of the source is lower (higher) than an energy of a relevant 
quantum state in the dot and, therefore, electrons tunnel 
from(to) the quantum dot through the barrier connected 
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